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Abstract 

State  feedback  regulators  are  derived  which  have  better  robustness 
characteristics  than  the  standard  60°  phase  margin  and  50%  gain  reduction 
tolerance  of  standard  linear-  quadratic  regulators.  It  is  also  shown  how 
the  Lyapunov  equation  can  be  used  to  design  high- integrity  regulators  for 


open- loop  stable  systems. 


1 .  Introduction 


Any  practical  control  system  synthesis  is  subject  to  uncertainty.  This 
uncertainty  may  appear  as  non-measurable  disturbances  on  the  input,  or  only 
partly  known  or  time  varying  system  parameters , to  quote  a  few  examples.  The 
topic  of  the  present  paper  is  to  study  one  example  of  such  parameter  un¬ 
certainty,  namely  large  variations  in  the  input  channels. 

The  criterion  chosen  for  our  "acceptable"  closed-loop  system  behavior 
is  stability.  This  is,  of  course,  far  from  sufficient  for  practical  purposes, 
but  such  analysis  may  all  the  same  produce  some  guidelines  for  the  design  of 
practical  control  sytems.  The  step  from  a  qualitative  to  a  quantitative  result 
is  often  a  short  one. 

Since  the  discovery  by  Kalman  [1]  of  the  frequency  domain  inequality 
satisfied  by  linear-quadratic  optimal  regulators,  a  great  deal  of  literature 
has  been  produced  on  the  "robustness"  of  such  regulators.  In  terms  of  classical 
control  concepts,  Kalman's  inequality  implies  that  they  possess  60°  phase 
margin,  infinite  gain  margin,  and  50%  gain  reduction  tolerance.  This  result 
holds  under  the  only  assumption  that  the  penalty  matrix  of  the  state  variable 
in  the  performance  index  is  positive  semide  initive. 

Although  60°  phase  margin  seems  a  lot,  it  may  be  insufficient  in  certain 
cases.  Dynamics  of  the  actuators  x->c  and  time  delays  may  have  been  neglected  in 
the  model,  for  instance.  Thus  it  is  of  interest  to  synthesize  regulators  with 
improved  robustness  characteristics. 

As  is  easy  to  see  from  examples,  the  above  results  cannot  be  strengthened 


without  further  assumptions  on  the  performance  index  or  the  plant.  Assumptions 
used  here  are  decreasing  penalty  on  the  control  variables  (  cheap  control) 
and  stability  of  the  open  loop  system,  respectively. 
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2.  Formulation  of  the  Problem. 

Consider  the  linear  time  invariant  system  described  by 
x(t)  =  Ax (t)  +  Bu (t) 

x  (0)  *  x  .  (1) 

o 

x (t)  and  u(t)  are  n-  and  p-vectors,  respectively.  Assume  that  the  desired  input 

is  given  by  the  linear,  constant  state  feedback 

u (t)  =  -  LTx(t)  (2) 

but  that  the  implementation  is  computed  by  some  time  varying  nonlinearity 

i.e.  the  actual  input  u(t)  is  given  by 

u(t)  -  -<J)(LTx(t)  ,t)  (2)'. 

P  P 

Here,  <N*»t)  is  a  nonlinear  function  from  R  to  R  for  each  t,  subject  to 
the  condition 

KjOTcj  <.aT<|>  fo,t)  <  VGe  RP  (3). 

The  problem  is  to  guarantee  global  asymptotic  stability  of  the  null  solution 
of  (1),  (2)'  subject  to  (3)  for  given  ,  K2- 

Although  the  problem  is  formulated  for  the  case  of  a  memoryless  non¬ 
linearity  4>(*,t) ,  the  results  can  easilty  be  generalized  to  dynamic  disturbances 
using  modern  frequency  domain  stability  theorems  (see  e.g.  [2]),  This  has  been 
exploited  in  [3] ,  but  will  not  be  pursued  here. 


For  notational  simplicity,  R  =  1^  will  be  assumed  in  the  sequel. 

Standard  manipulations  of  equation  (5) lead  to  the  frequency  domain 
equality 

(I  +  GT(-ilo))I  (I  +  G(iu>) )  =  I  +  BT(-idl  -  aV1  Q  (iw  -  A)"^  (6), 

P  P 

where 

G(s)  =  R-1BTP(SI  -  A)_1B. 

An  alternative  form  of  (6)  for  the  closed-loop  system  is,  with  G(s)  -  G(s)(I+G(s)) 
(I  -  GT(-iu)))Ip(l  -  G(iu) )  =  I  -  BT(-iuj  -  AT)_1  Q(iu>  -  A)_1B  (6)'. 

If  Q  is  positive  sem:  definite ,  the  second  term  on  the  right-hand  side  is 
nonnegative,  and  the  following  inequality  results: 

(I  -  GT(-iu>))I  (I  -  G (i<D) )  <  I  (7). 

r  r 

A  straight-forward  application  of  the  circle  criterion  then  implies  the  desired 
stability  for 

Ki  =  j,  =  any  finite  K. 

It  is  obvious  that  the  best  possible  lower  bound  without  restrictions  on  the 
plant  is 

Kx  >  0  (8). 

For  certain  high-gain  regulators,  (8)  is  in  fact  sufficient  as  will  now  be 
shown.  Consider  the  performance  index  Jp,  where  the  penalty  on  the  input  is 
modified  via  the  scalar  variable  p: 

Jp  -  J  (x(t)T  Qx(t)  +  pu(t)T  u(t))dt. 

Choosing  p  small  means  that  very  large  control  signals  are  permitted.  It 
seems  reasonable  to  conjecture  that  such  a  design  should  yield  a  controller 
with  high  integrity.  This  is  indeed  the  case. 

Theorem  1.  Consider  the  control  system  (1),  (2),  with  ♦(•,t)  subject  to  (3). 


For  any  >0  and  any  K2  <  00  there  is  a  Q  ii  0  such  that  the  optimal  controller 
relative  to  Jp  quarantees  asymptotic  stability  in  the  large  of  the  null  solution 
for  all  p  >0  sufficiently  small. 

Proofs  The  Riccati  equation  pertaining  to  is 

ATPp  +  PpA  +  Q  -  P_1PpB  BTPp  =  0  (5)  '  . 

Two  reductions  of  the  general  case  will  be  made.  Firstly ,  it  may  assumed  without 
loss  of  generality  that 

rank  (Q)  £  rank  (B)  =  p, 

T 

(see  [4]).  Thus  Q  can  be  factored  as  MM  ,  with 
rank  (M)  <  rank  (B) . 

There  are  several  M-matrices  that  generate  the  same  feedback  matrix, 

corresponding  to  various  factorizations  of 

T  T  -1  -It  t  “1  t  -1 

B  (-si  -  A  )  Q(sl  -  A)  B  =  B  (-si  -  A  )  MM  (si  -  A)  B. 

If  the  M  that  makes 

M(s)  =  MT(SI  -  A)"1  B 

minimum  phase  is  chosen,  the  conditions  are  satisfied  which  ensure  that  (15], 

[6]) 

lim  P  *  0 
P-K>  P 

Equivalently, 

lim  p_1PpB  BT  P  =  Q  =  MMT  (9) . 

p**o 

Starting  from  Riccati  equation  (5) 1  and  using  the  same  manipulations  that 
lead  to  inequality  (6)’  then  shows,  inserting  (9),  that 

G  (s)  *  p“1BTPp(sl  -  A  +  P_;LBBTPp)“1  B 
satisfies 

<jl  -  GpT(-s))Ip(il  -  Gp(s))  <Ip  +  0(p_1),  p  +  0. 

This  proves  the  claim. 


Remark:  A  weaker  version  of  Theorem  1,  corresponding  to  a  4>(*,t)  that  is 


diagonal  (i.e.  no  cross-couplings  between  the  inputs),  can  be  proved  using  the 
results  of  [7] . 


In  classical  terminology,  this  regulator  possesses  a  phase  margin  arbitrary 
close  to  90°,  infinite  gain  margin,  and  a  gain  reduction  tolerance  arbitrarily 
close  to  100%. 

There  is  an  alternative  way  of  generating  such  controllers,  which  is 
almost  trivial.  The  proof  is  left  to  the  reader. 

Theorem  2.  Consider  (1),  (2)'  with  <t>  (•  ,t)  subject  to  (3).  If  L  is  chosen  as 

L  =  G*1  L  , 
o 

where  Lq  is  any  optimal  gain  generated  from  a  positive  semidefinite  Q,  then 
the  controller  enjoys  the  same  robustness  properties  when  G  +  0  as  the  controller 
of  Theorem  1  when  p  0. 

Theorems  1  and  2  must  be  used  with  some  caution,  since  in  applications  the 
sector  bounds  and  K2  may  depend  on  the  nominal  feedback  gain  L.  An  example  is 
given  in  Section  5. 

4.  Results  for  Open-Loop  Stable  Plank 

If  the  plant  is  open-loop  stable,  it  should  be  possible  to  design  controllers 
where  the  lower  bound  is  zero.  A  complete  characterization  of  such  feedback 
matrices  L  is  given  in  the  following  theorem. 

Theorem  3.  The  feedback  system  (1),  (2)'  with  4>(*, t)  satisfying  (3)  has  a 
globally  asymptotically  stable  equilibrium  solution  for 
Ki  =  0,  K2  =  any  K  >  0 

T 

if  there  exist  matrices  K  =  K  >0  and  C,  (C,A)  being  an  observable  pair,  such 
that 

/  T  T 

1 A  K  +  KA  =  -CC 


KB  *  L. 


(10) 


,'J mitx> 
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Proof:  This  is  a  straight-forward  application  of  the  Kalman  Yakubovich- 

Popov  lemma  in  its  multivariable  form  ([8]).  Equation  (10)  provides  a  means 
to  design  high- integrity  controllers,  namely  by  solving  the  (linear)  Lyapunov 
equation  for  K  and  then  choosing  the  feedback  gain  L  as  (any  positive  multiple  of) 
KB.  There  remains  the  problem  of  choosing  C  and  a  positive  gain,  but  this 
problem  is  shared  by  the  LQOC  methodology.  For  single-input  systems,  equation  (10) 
is  simplified  further  if  it  is  translated  into  the  frequency  domain.  Defining 
p(s)  as  the  open-loop  characteristic  polynomial  and 


q(s) 

P(s) 


T  -1 

L  (si  -  A)  B 


yields  the  equation 


A)”1  B 


p(s)  q(-s)  +  q(s)  p(-s)  =  r(s)  r(-s)  (10)’  . 

This  can  be  solved  for  q(s),  which  is  turn  gives  the  unique  feedback  gain  L. 

The  condition  of  Theorem  3  is  rather  strong.  In  general  an  infinite 
gain  margin  is  not  required.  The  following  theorem  characterizes  the  feed¬ 
back  gains  that  retain  stability  for  all  gain  drops. 


Theorem  4.  The  claim  of  Theorem  3  remains  valid  for 
K1  *  0r  *2  ~  1 

T 

if  there  exist  matrices  K  *  K  >  0,  C,  and  D,  and  an  e  >  0  such  that 

T  T  T 

A  K  +  KA  «  -  CC  -  DD 

(id 

KB  *  1  -  (/2  -  e)D 


Proof.  This  is  again  a  direct  application  of  the  circle  criterion  in  combination 
with  the  Kalman- Yakubovich-Popov  lemma.  Compared  to  equation  (10),  the  use  of 
equation  (11)  involves  the  extra  problem  of  choosing  D.  Theorem  4  can  of  course 
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be  generalized  to  any  upper  and  lower  bounds  and  (provided  the  system 
is  open  loop  asymptotically  stable  if  <_  0  <_  K2- 

5.  Discussion  of  the  Results 

The  results  presented  above  must  be  interpreted  with  some  caution.  It  is 
quite  clear  that  most  actuator  failures  are  such  that  no  judicious  choice  of 
feedback  gains  can  save  the  situation.  This  is  even  more  true  for  the  case 
of  sensor  failures  (not  treated  here) .  Thus  in  general  one  will  have  to  rely 
on  external  failure-handling  routines. 

The  sharpening  of  previous  results  for  arbitrary  plants  relies  upon 
the  use  of  high  gains.  This  must  be  kept  in  mind  if  the  nominal  input  influences 
the  sector  of  4>(*,t).  For  instance,  if  the  gain  drop  of  the  regulator  is  caused 
by  saturation,  there  is  obviously  no  point  in  increasing  the  nominal  gain. 

With  these  reservations  kept  in  mind,  the  theorems  yet  contribute  to 
the  general  robustness  picture  of  state  feedback  regulators.  Theorem  3  is 
believed  to  yield  a  practicable  design  method  for  control  systems  where  the 
input  channel  uncertainty  is  so  great  that  it  has  to  be  accounted  for  in  the 


synthesis. 
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